Origin of peraluminous granites and granodiorites, Iberian massif, Spain: an experimental test of granite petrogenesis

The discrimination between potential source materials involved in the genesis of Iberian granites and granodiorites, as well as the role of mantle-crust interactions, are examined using constraints imposed by melting experiments, melting-assimilation experiments and Sr-Nd isotope systematics. The Sr-Nd isotope relationships indicate the existence of different genetic trends in which juvenile mantle materials are involved by different mechanisms: (1) a source trend, traced by a particular evolution of the pre-Hercynian basement and indicating mantle participation at the time of sedimentation; (2) a set of magmatic trends traced by gabbro-tonalite-enclave-granodiorite associations, implying the incorporation of new mantle material at the time of granite generation. These relationships strongly support a pure crustal origin for the peraluminous leucogranites, derived from partial melting of crustal protoliths, and a hybrid origin for the peraluminous granodiorites. These granodiorites are the most abundant granitic rocks of the Central Iberian zone (CIZ) of the Iberian massif, implying that processes of hybridisation by assimilation and/or magma mixing played an important role in granitoid production during the Hercynian orogeny. These hypotheses have been tested by means of melting and assimilation experiments. Melting experiments in the range 800–900 °C and at pressures of 3, 6, 10 and 15 kbar indicate that: (1) several potential source materials such as Bt-Ms gneisses and metagreywackes are suitable for the production of peraluminous leucogranite melts; (2) the melt compositions are always leucogranitic, regardless of pressure; (3) pressure exerts a strong influence on the fertility of the source: experiments at 3 kbar produce more than 20 vol% of melt, compared with less than 5 vol% of melt produced at 10 and 15 kbar and at the same temperature. The melting-assimilation experiments carried out at 1000 °C and 4, 7 and 10 kbar and using a proportion of 50% gabbro and 50% gneiss give high melt proportions (more than 50 vol.%) and noritic residues. These melts have the composition of leucogranodiorites, and overlap with part of the compositional range of peraluminous granodiorites of the Iberian massif. The generation of more mafic granodiorites may be explained by the incorporation of some residual orthopyroxene to the granodiorite magmas. The low solubility of Fe + Mg prevents the generation of granodiorite melts with more than 3 wt% of MgO + FeO at all crustal pressures. The large volumes of peraluminous, hybrid granodiorites, produced by assimilation of crustal rocks by mantle magmas, imply that an important episode of crustal growth took place during the Late-Palaeozoic Hercynian orogeny in the Iberian massif. Received: 30 June 1998 / Accepted: 27 November 1998     

The effect of seafloor topography in the Southern Ocean on tabular iceberg drifting and grounding

Antarctic tabular icebergs are important active components in the ice sheet-ice shelf-ocean system. Seafloor topography is the key factor that affects the drifting and grounding of icebergs, but it has not been fully investigated. This study analyzes the impact of seafloor topography on the drifting and grounding of Antarctic tabular icebergs using Bedmap-2 datasets and iceberg route tracking data from Brigham Young University. The results highlight the following points. (1) The quantitative distributions of iceberg grounding events and the tracking points of grounded icebergs are mainly affected by iceberg draft and reach their peak values in sea water with depths between 200 m and 300 m. The peak tracking point number and linear velocity of free-drifting icebergs are found in the Antarctic Slope Front (water depth of approximately 500 m). (2) The area of possible grounding regions of small-scale icebergs calved from ice shelf fronts accounts for 28% of the sea area at water depths less than 2000 m outside the Antarctic coastline periphery (3.62 million km²). Their spatial distribution is mainly around East Antarctica and the Antarctic Peninsula. The area of possible grounding regions of large tabular icebergs with long axes larger than 18.5 km (in water depths of less than 800 m) accounts for 74% of the sea area. (3) The iceberg drifting velocity is positively correlated with ocean depth in areas where the depth is less than 2000 m (R=0.85, P<0.01). This result confirms the effect of water depth variations induced by seafloor topography fluctuations on iceberg drifting velocity.     

Enhancing the Forecasting of Monthly Streamflow in the Main Key Stations of the River Nile Basin

Predicting the streamflow of rivers can have a significant economic impact, as this can help in agricultural water management and in providing protection from water shortages and possible flood damage. In this study, two statistical models have been used; Deseasonalized Autoregressive moving average model (DARMA) and Artificial Neural Network (ANN) to predict monthly streamflow which important for reservoir operation policy using different time scale, monthly and 1/3 monthly (ten-days) flow data for River Nile basin at five key stations. The streamflow series is deseasonalized at different time scale and then an appropriate nonseasonal stochastic DARMA (p, q) models are built by using the plots of Partial Auto Correlation Function (PACF) to determine the order (p) of DARMA model. Then the deseasonalized data for key stations are used as input to ANN models with lags equals to the order (p) of DARMA model. The performance of ANN and DARMA models are compared using statistical methods. The results show that the developed model (using 1/3 monthly (ten-days) and ANN) has the best performance to predict monthly streamflow at all key stations. The results also show that the relative error in the developed model result did not exceed 9% while in the traditional models reach to 68% in the flood months in the testing period. The result also indicates that ANN has considerable potential for river flow forecasting.     

Assessment of models predicting anthropogenic interventions and climate variability on surface runoff of the Lower Zab River

Multi-regression, hydrologic sensitivity and hydrologic model simulations were applied to quantify the climate change and anthropogenic intervention impacts on the Lower Zab River basin (LZRB). The Pettitt, precipitation-runoff double cumulative curve (PR-DCC) and Mann–Kendall methods were used for the change points and significant trend analyses in the annual streamflow. The long-term runoff series from 1979 to 2013 was first divided into two main periods: a baseline (1979–1997) and an anthropogenic intervention period (1998–2013). The findings show that the mean annual streamflow changes were consistent using the three methods. In addition, climate variability was the main driver, which led to streamflow reduction with contributions of 66–97% during 2003–2013, whereas anthropogenic interventions caused reductions of 4–34%. Moreover, to enhance the multi-model combination concept and explore the simple average method (SAM), Hydrologiska Byrans Vattenbalansavdelning (HBV), Génie Rural a Daily 4 parameters (GR4J) and Medbasin models have been successfully applied.     

A Dual Delineation Approach to Untangle Interferences in Well Doublets Systems

Traditional numerical methods for the delineation of wellhead protection areas span deterministic and probabilistic approaches. They provide time-related capture zones. However, none of the existing approaches identifies the groundwater contribution areas related to each source or sink. In this work, the worthiness of the so-called double delineation approach was extended. This task was achieved by simple postprocessing of its dual outputs leading to a highly efficient screening tool. In the particular context of geothermal resources management through the well doublets of the Dogger aquifer in the Paris Basin (France), the approach was extended to forecast the compositional heat breakthrough at production wells. Hence, cold-water breakthrough and temperature decline in production wells are timely assessed in low-enthalpy geothermal reservoirs. The method quantifies how groundwater volumes are moving through space and time between any couple of source and sink. It provides unprecedented tools advancing the enhanced understanding of water resources systems functioning. It is highly recommended to implement the presented concepts in the current and future generations of community groundwater models.     

Lithology-dependent seismic anisotropic amplitude variation with offset correction in transversely isotropic media

Analysis of amplitude variation with offset is an essential step for reservoir characterization. For an accurate reservoir characterization, the amplitude obtained with an isotropic assumption of the reservoir must be corrected for the anisotropic effects. The objective is seismic anisotropic amplitude correction in an effective medium, and, to this end, values and signs of anisotropic parameter differences (Δδ and Δε) across the reflection interfaces are needed. These parameters can be identified by seismic and well log data. A new technique for anisotropic amplitude correction was developed to modify amplitude changes in seismic data in transversely isotropic media with a vertical axis of symmetry. The results show that characteristics of pre-stack seismic data, that is, amplitude variation with offset gradient, can be potentially related to the sign of anisotropic parameter differences (Δδ and Δε) between two layers of the reflection boundary. The proposed methodology is designed to attain a proper fit between modelled and observed amplitude variation with offset responses, after anisotropic correction, for all possible lithofacies at the reservoir boundary. We first estimate anisotropic parameters, that is, δ and ε, away from the wells through Backus averaging of elastic properties resulted from the first pass of isotropic pre-stack seismic inversion, on input data with no amplitude correction. Next, we estimate the anisotropic parameter differences at reflection interfaces (values and signs of Δδ and Δε). We then generate seismic angle gather data after anisotropic amplitude correction using Rüger's equation for the P-P reflection coefficient. The second pass of isotropic pre-stack seismic inversion is then performed on the amplitude-corrected data, and elastic properties are estimated. Final outcome demonstrates how introduced methodology helps to reduce the uncertainty of elastic property prediction. Pre-stack seismic inversion on amplitude-corrected seismic data results in more accurate elastic property prediction than what can be obtained from non-corrected data. Moreover, a new anisotropy attribute (ν) is presented for improvement of lithology identification.     

Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Research note: deblended-data reconstruction using generalized blending and deblending models

We introduce a concept of generalized blending and deblending, develop its models and accordingly establish a method of deblended-data reconstruction using these models. The generalized models can handle real situations by including random encoding into the generalized operators both in the space and time domain, and both at the source and receiver side. We consider an iterative optimization scheme using a closed-loop approach with the generalized blending and deblending models, in which the former works for the forward modelling and the latter for the inverse modelling in the closed loop. We applied our method to existing real data acquired in Abu Dhabi. The results show that our method succeeded to fully reconstruct deblended data even from the fully generalized, thus quite complicated blended data. We discuss the complexity of blending properties on the deblending performance. In addition, we discuss the applicability to time-lapse seismic monitoring as it ensures high repeatability of the surveys. Conclusively, we should acquire blended data and reconstruct deblended data without serious problems but with the benefit of blended acquisition.     

Stress distribution from railway track over geogrid reinforced ballast underlain by clay

Railway ballast forms a major component of a conventional rail track and is used to distribute the load to the subgrade, providing a smooth running surface for trains. It plays a significant role in providing support for the rail track base and distributing the load to the weaker layer underneath. Ballast also helps with drainage, which is an important factor for any type of transportation structure, including railroads. Over time, ballast progressively deforms and degrades under dynamic loading and loses its strength. In this study, extensive laboratory tests were conducted to investigate the effect of load amplitude, geogrid position, and number of geogrid layers, thickness of ballast layer and clay stiffness on the behavior of the reinforced ballast layer and induced strains in a geogrid. A half full-scale railway was constructed for carrying out the tests, which consisted of two rails 800 mm in length with three wooden sleepers(900 mm × 10 mm × 10 mm). Three ballast thicknesses of 200, 300 and 400 mm were used in the tests. The ballast was overlying 500 mm thickness clay in two states, soft and stiff. The tests were carried out with and without geogrid reinforcement; the tests were performed in a well-tied steel box of 1.5 m length ×1 m width ×1 m height. Laboratory tests were conducted to investigate the response of the ballast and the clay layers where the ballast was reinforced by a geogrid. Settlement in ballast and clay, soil pressure and pore water pressure induced in the clay were measured in reinforced and unreinforced ballast cases. It was concluded that the amount of settlement increased as the simulated train load amplitude increased, and there was a sharp increase in settlement up to cycle 500. After that, there was a gradual increase that leveled out between, 2500 to 4500 cycles depending on the frequency used. There was a slight increase in the induced settlement when the load amplitude increased from 0.5 to 1 ton but it was higher when the load amplitude increased to 2 tons. The increased amount in settlement depended on the existence of the geogrid and other parameters studied. The transmitted average vertical stress for ballast thicknesses of 30 cm and 40 cm increased as the load amplitude increased, regardless of the ballast reinforcement for both soft and stiff clay. The position of the geogrid had no significant effect on the transmitted stresses. The value of the soil pressure and pore water pressure on ballast thicknesses of 20 cm was higher than for 30 cm and 40 cm thicknesses. This meant that the ballast attenuated the induced waves. The soil pressure and pore water pressure for reinforced and unreinforced ballast was higher in stiff clay than in soft clay.